The violation of the Hund's rule in semiconductor artificial atoms 
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The unrestricted Pople-Nesbet approach for real atoms is adapted to quantum dots, the man- 
made artificial atoms, under applied magnetic field. Gaussian basis sets are used instead of the exact 
single-particle orbitals in the construction of the appropriated Slater determinants. Both system 
chemical potential and charging energy are calculated, as also the expected values for total and 
z-component in spin states. We have verified the validity of the energy shell structure as well as the 
Hund's rule state population at zero magnetic field. Above given fields, we have observed a violation 
of the Hund's rule by the suppression of triplets and quartets states at the lp energy shell, taken as 
an example. We also compare our present results with those obtained using the LS-coupling scheme 
for low electronic occupations. We have focused our attention to ground-state properties for GaAs 
quantum dots populated up to forty electrons. 
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The influence of spatial confinement on the physical 
properties such as the electronic spectra of 

0D structures is a topic of growing interest. Among 
such systems one could select carriers and impurity atoms 
in metallic or semiconductor mesoscopic structures^ as 
also atoms, ions and molecules trapped to microscopic 
cavities (i2i2i^. In these systems, the confinement be- 
comes important whenever a quantum sizes equal the 
cavity length. However, the energy spectrum of these 
systems is not only determined by the spatial confine- 
ment and its geometrical shape, but also by environmen- 
tal facts such as electric and magnetic applied fields that 
break or lower the general symmetries. Finally, many- 
body effects as electron-electron interaction, may be even 
be more important than the confinement itself. In any 
case, a correct description of physical properties of the 
problem requires that the wavefunction must reflect both 
the form of confinement and the appropriated boundary 
conditions. 

Interesting confined systems are the semiconductor 
quantum dots structures (QDs), also referred as man- 
made artificial atoms, built as low-dimensional electronic 
gases when crystalline translation invariance is broken 
in all three spatial dimensions and leading to discrete 
energy states, as in real atoms. Various are the ap- 
proaches that have been used to deal with many-particle 
QDs. Among them, one can cite charging model, 6 ' 7,8,9 
correlated electron model^ Green functions^ Lanczos 
algorithm^ Monte Carlo method^ Hartree-Fock (HF) 
calculations fi^i^iSiii and density functional theory^ 

The charging model can reproduce well experimental 
findings for "metallic dots". On the opposite side, with 
much lower electronic density, the semiconductor dots 
requires a microscopic point of view to treat the electron- 
electron interaction. 

Here we will consider a QD defined by a "hard wall" 
spherical volume which is appropriated for semiconduc- 



tors grown inside glass matrices. Some of the commonly 
studied topics are the formation of energy shells in their 



spectra^ the control of their electronic correlations 
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the formation of Wigner moleculeSfSi and the influence 
of the Coulomb interaction in their spectrai 2 ^ 3 . In these 
spherically defined artificial atoms both spin and orbital 
angular momenta are good quantum numbers, and the 
low occupation many-particle eigenstates can be labelled 
according to the LS-coupling scheme? 2 ^. For occupation 
number above four, the LS-coupling scheme becomes ex- 
tremely cumbersome and, in this paper, we have chosen 
the unrestricted Pople-Nesbet matrix approach 2 ^ of the 
single determinant self-consistent HF formalism to treat 
shell configurations of dots containing up to forty elec- 
trons, where we have calculated the total spin expected 
values, chemical potential and charging energy. We show 
the changes induced by the magnetic field on such ap- 
proach, using a set of Gaussian basis (section II). Then 
we discuss our main results (section III) by focusing on 
how magnetic field determines the Zeeman splitting and 
induces violation of the Hund's rule. 

Within the Unrestricted Hartree-Fock formalism 
(UHF), the a (spin- up) and f3 (spin-down) func- 
tions have different spatial components, Xi( x ) = 

described by the orbitals 
1. ..&}). Therefore, an UHF 



{^ Q (r)a( W ),^(r)/3(a;)}, 

{'.'/./ i •••/••!• = 



i a i P 



which 



wavefunction has the form \ty UHF } 

represents open shells once no spatial orbital can be dou- 
bly occupied. The closed shells can be also obtained, 25 
more specifically, UHF functions are not necessarily sys- 
tem eigenstates having well defined L and S values. Yet, 
the number of carrier, N , must equal the the sum of spin- 
up and spin-down electrons, as N = N a + N@ . The inte- 
gration of the spin degrees of freedom 2 ^ yields two cou- 
pled HF equations that must be simultaneously solved. 
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They have the form v - j 

the respective Fock operators are given by 



a/B\ 



a/3 



ip°[ ), where 



P, 



IP 



/J W 



/B^a/3* 



Thus, 



one can 



de- 



fine two new quantities: i) The total charge density, 
p T (r) = p a (r) + p@{r), that determines N when in- 



F«x/P 



tegrated over all space; ii) The spin density, p (r) 



- K2 /l3 \ + E J% /a . (f) P"( r ) ~ P fi ( v )i tnat yields 2M S after integration over all 
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Both f a and f 13 include the kinetic (hj), the direct( J" ) 
and the exchange (K"^) terms between electrons with 

same spin, and a direct term (Ja^ a ) between electrons 
with opposite spins. The interdependence among f a (f@) 
and 

requires the simultaneous solution of the two 

HF equations. They yield the sets {1/1™ } and (ij>^ j that 

should minimize the energy Eq HF of the unrestricted 
ground-state, |\f , o' ffF ), given by 



E 



UHF 



N" 



1 I'D/I 



JV' 3 N a N" 

E^ + oEE^-^r] 
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N f3 N a N 13 
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In these expressions, hai^ 



/0 



ha 



I a/8 



J 



aB 
ab 



= J, 



pa 



jaa 

J ab — 



(CI^IC) = {< Ja 

(CI 4* IC> = (CI 4" IC> = J feT (with an analog term 
for J%), and ^ = (CI IC> = (CI^IC) = 
-^6a* ( w ith an analog term for K^). 

The Pople-Nesbet approach transforms the UHF equa- 
tions into a matrix formulation by expanding ipf^ in a 
set of known basis functions {4> u \i> = !...&}, 



space. The unrestricted wavefunctions are eigenfunctions 
of Sz-, but not necessarily of S 2 , therefore one can de- 
fine the total charge ( P T = P" + P' 3 ) and the spin ( 
pS _ pa _ p/3-j d ens ity matrices for the system. The el- 
ements of the two Fock matrices are obtained as F"/^ = 



a/3 



V" \<t> v ) and 

\4> v <t>\) - 



e //, where T M „ = -h 2 /(2m) 

\f (<f>»<t>a\ l r i - r 2| _1 I^a<^)], with s being the mate- 
rial dielectric constant. 

The self-consistency lies in the fact that both F and P 
depend on C, while the coupling of spin-up and -down 
equations occurs since F a (F' 3 ) depends on P@ (P Q ) 
through P T . 

The procedure to solve Eq. 0]is: i) Given a confinement 
potential, one specifies ./V and {^j; ii) The integrations 
on S^v and T^ v are performed; iii) An initial guess is 
used for P"^ and P T , perform two-electron integrals for 
G a/(3 and T to construct F a ^ . Diagonalize F"/" to get 
Qa/g an( j £ a/f) ^ an( j f orm new pa//3- j v ) This iteration 

is repeated until the desired convergence for E^ HF is 
reached. The Pople-Nesbet ground-state energy is 



j^UHF 



9 E E 



pa pa , p p3 1 /gN 



The unrestricted functions are not, in general, eigen- 
states of S 2 , thus we calculate the spin expectation values 



.,25 



a /g _ 



E^; 



(3) 



(S 2 ) 



UHF 



N a 



N 



where the expansion coefficients become the pa- 

rameters to be iterated. When Eq. [31 is inserted into Eq. 
^one obtains the k x k coupled matrix equations, 



pa//3(-iQ//3 _ SQ a / P £ a /P 



(4) 



where S is the positive defined overlap (5*^^ = (<^>„ \(j> v )) 
between basis functions, C"/' 9 are the expansion coeffi- 
cient matrices whose columns describe each spatial or- 
bital are the diagonal matrices of the orbital 
energies (e™ ), and F"/* 3 are the Fock operators, 



F^ = / «(r)f^„(r). 



(5) 



At this point becomes convenient to introduce the charge 
density for the spin-up and -down electrons, defined as 
p a/ 3(r) = E JV«"> | C /^ (r) | 2 = E.E.^^WCW, 
where the elements of the respective density matrices 



N N 

-EE 



and 



(Sz) 



UHF 



E E (■^ • 



(7) 



(8) 



As an application of the Pople-Nesbet approach, we 
consider a QD with radius Rq confined to an infinite 
spherical potential in the presence of a magnetic field 
B = (0,0, Bo) and populated up to forty electrons. Its 
single-particle Hamiltonian is 



H 



2m 



^ 2 
lie 



S ^B.S, 



(9) 



where p B is the Bohr magneton, g is the bulk g-factor, 
and we use the symmetric gauge, A = (B x r)/2. Using 
atomic units, Er v — e 2 /(2a )(Rydberg) for the energy, 



3 



and do = /i 2 /(moe 2 )(Bohr radius) for length, the Hamil- 
tonian H can be written in dimensionless form 



l_o| 
mR 2 



x 2 dx 



,d_ 

dx 



L 2 

x 2 



(10) 



R 2 
21 



mgS 2 



where rh = m/mo, Ib — y/hc/(eBo) is the magnetic 
length, and x = r/Ro. Without magnetic field, the nor- 
malized spatial eigenfunctions of Hq are given by 



0» 



1 



[ji+i (a n i)Y 



1/2 



jl (oinix)Yi <mi (6,i 
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The boundary condition at the surface, r — Rq (or 
x = 1), determines a n i as the n th zero of the spherical 
Bessel function ji(a n ix) and Yi tmi (8,(f>) is the spherical 
harmonic. 

The Hamiltonian for the electron-electron interaction 
in atomic units becomes H ee = (ao/Ro)2/(e |xi — X2I), 
where the usual multipole expansion for |xi — X2I -1 is 
used in our calculations. 

The spatial orbitals in ipt^ define the six lowest energy 
shells (Is, lp, Id, 2s, If, 2p) without magnetic field^. 
Thus, the index v = n,l,mi can assume up to 40 (20 
spin-up and 20 -down) possible values for those shells. 
Certainly, the magnetic field lifts both spin and orbital 
degeneracies. Let us consider a GaAs dot, a wide-gap 
semiconductor having rh — 0.065, g = 0.45 and e = 
12.65. 

The inclusion of a magnetic field requires modifications 
on the UHF equations. The spin-independent linear and 
quadratic magnetic terms are easyly added to the defini- 
tions of both h in f a '@ and T Mt/ . However, the inclusion of 
the spin-dependent linear term (~ BqSz) to h in f a /P re- 
quires decomposition of kinetic and Coulomb terms into 
'"" 1 Thus, under a magnetic field we should 

in 



Tfj,u H and G^ v 

make the substitution P^T^ 



i/fi [if 



P/3 rpfi 



E% HF (Eq.EJ. 

Another important detail refers to the orbital basis 
{(j) v \v = l...fc}. Instead of the exact spherical Bessel func- 
tions of Ea llll the radial part of each orbital was decom- 
posed in a sum involving five Gaussians confined to region 
x < 1 , while the angular part is maintained as defined by 
its symmetry, as 

n-i 

Ki,m t (x,e,4>) = Nm(l-x) n x l l[(au-x) (12) 

1=1 
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- DkR ° x Y l>mi {6,cj>), 



where N n i is the normalization, the polynomial (1 — x) 
garanties the boundary x = 1, the polynomial in x l is 
required for I > states at the origin x — 0, the product 



? 20 _ GaAs - R = 1 00 A - B fl = 
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% — o — Non-interacting Electrons 
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FIG. 1: Unrestricted ground-state energy for a Ro = 100 A 
GaAs QD without magnetic field. In the upper left panel we 
compare the unrestricted and the non-interacting energies, 
where the QD energy shell structure is visible. The bottom 
left panel shows QD chemical potential (left scale) and charg- 
ing energy (right scale). The former displays abrupt change 
always when a new shell starts to be populated, while the 
latter presents larger (smaller) peaks when a shell is fully 
filled (half-filled), a direct consequence of the Hund's rule. 
The bottom and upper right panels show, respectively, the 
iV-evolution of the expected values of total spin and its z- 
projection. 



in (an — x) nulls functions at the zeros an of the respec- 
tive spherical Bessel function transposed to the interval 
< x < 1, and the last sum involves an expansion into 5 
Gaussians. Higher order expansion did not show any im- 
provement for N < 40. The Gaussian coefficients 14 and 
exponents are determined for each value of Rq, by 
maximizing the superposition between Eqs. El an d ITSl 
Once 14 and Dk are determined, we run our UHF code 
for each value of i?o and N , and find the parameters 
C"/' 3 that better describe Eq. [3] and give the minimal 
energy in Eq. |SJ 

At last, we have calculated two quantities that will be 
used in the description of our results. The first one is the 
QD chemical potential, which yields the energy difference 
between two successive ground states, 



tx dot (N) = E (N)-E (N-l) 



(13) 



The second one is the QD charging energy, which yields 
the energy cost to add an extra electron to the system, 



E 



(N) = E (N + 1) — 2E (N) + E (N — 1) . (14) 



'char 



From these two last equations, one can also see that 
E char (N) = ii dot (N +1) - fi dot (N). 

We show in Fig. ^the results of a UHF Pople-Nesbet 
calculation for a GaAs QD having Rq — 100 A, at zero 
magnetic field. In the left upper panel we have compared 
the non-interacting electron problem and the UHF results 
as function of occupation. The shell structure occurs 
for magic numbers N = 2, 8, 18, 20, 34 and 40. It 
is observed that electron-electron interaction decreases 
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FIG. 2: Violation of Hund's rule induced by magnetic field in 
the Ro = 100 A GaAs QD of the previous figure. The panels 
show the successive occupation (indicated in the upper right 
corner of each panel) of the lp shell, assuming that the Is 
shell remains populated by one spin-up and one spin-down 
electron. The possible spin configurations for given N are 
indicated by + (spin-up) and — (spin-down). For Bo = the 
spin sequence is 1/2 — 1 — 3/2 — 1 — 1/2 — 0, while at fields 
higher than 3 T, it changes to 1/2 - - 1/2 - - 1/2 - 0. It 
interesting to observe that a magnetic field can suppress the 
energy triplets and quartets from the QD spectrum. 

(increases) the non-interacting ground-state energy when 
occupation corresponds to a shell less (more) than half- 
filled. At exactly half-filled cases N — 5, 13, 27 and 
37, the interacting and non-interacting energy values are 
approximately equal. 

In the bottom left panel of Fig. 2] we show both QD 
chemical potential (left scale, Eq. I13fl and charging en- 
ergy (right scale, Eq. I14fl . where the respective values 
of Eq were obtained from the unrestricted calculation 
presented in the left upper panel. Notice that n dot lin- 
early increases as the occupation increases inside a given 
shell. When such shell is totally filled, there is an abrupt 
change in fi dot indicating that the following shell starts 
its occupation. Also, the higher the occupation, the more 
abrupt is this change. An anomalous behavior seems to 
occur with the 2s shell, whose /i dot -value is larger than 
the 1/ shell, that has higher energy. The charging energy 
is another form to verify not only the presence of shell 
strusture in the spectrum, but also the validity of Hund's 
rule for the filling of such shells. In principle, E c h ar must 
present larger (smaller) peaks when the total (half) oc- 
cupation of a shell is achieved. The first fact is due to the 
higher difficulty in adding an electron to a QD in a filled 
shell state, while the second one refers to Hund's rule, 
which states that electrons must be added to the system 
with their spin being parallel, until all possible orbitals 
inside a given shell be occupied, making the total energy 
of the system be decreased because of the maximized ex- 
change contribution. However, some discrepancies are 
verified in E c h ar : the smaller peak of N = 27 occurs here 
at N = 26, and the larger peak of N — 20 is negative. 

The right bottom and upper panels of Fig^ show re- 



spectively the Af-evolution of the total spin S and its 
projection Ms as calculated from Eqs.[7|and|S|for the un- 
restricted energies. Notice that, with no magnetic field, 
Hund's rule seems to be followed; the Ms expected value 
oscillates from in a filled shell to its maximum in a half- 
filled shell, when it starts to decrease again on the way to 
the closing of the shell; the maxima are Ms = 1/2, 3/2, 
5/2 and 7/2 for s, p, d and / shells, respectively. The S 
expected value yielded by the unrestricted formalism is 
also very reasonable; discrepancies are only observed at 
iV = 24, where S > 2, and at N = 21, where S > 1/2. 
We believe that both discrepancies related to the 2s shell 
or to its surroundings - fi dot larger than the one of 1/ 
shell, negative peak for N = 20 in E c har, and almost 
doubled S expected value for N — 21 - are caused by the 
non-reasonable Gaussian reproduction of this orbital^i 

By focusing on the lp shell we show in Fig. |3 for the 
same QD of the previous figure, how a finite magnetic 
field is able to violate Hund's rule in the system. Pan- 
els from left to right and from up to bottom show the 
successive ground state energies from N = 3 to N = 8 
as this shell is filled, always considering that the Is shell 
remain fully occupied by two electrons, one up and one 
down; the distinct possible spin configurations for each 
N are indicated by + (up) and — (down). In addition to 
the small Zccman effect present in all occupations, there 
is a changing of ground state spins at N = 4, 5 and 6 
as the field is increased. Notice that at zero field the 
spin sequence is 1/2 — 1 — 3/2 — 1 — 1/2 — 0; in a field 
above 3 T it becomes 1/2 — — 1/2 — — 1/2 — 0, meaning 
that quartets and triplets are suppressed by the magnetic 
field, and the ground state of the system starts to oscil- 
late only between singlets and doublets at high fields as 
TV increases. When this lp shell is half- filled (N = 5), the 
ground state goes from a quartet to a doublet at Bq ~ 2 
T; when it has one electron more (N — 6) or less (N = 4) 
than that, it goes from a triplet to a singlet at Bq ~ 3 T. 

At last, in order to prove the efficiency of the Pople- 
Nesbet approach, we compared the results from this UHF 
self-consistent matrix formulation with the ones obtained 
from the LS-coupling scheme used in Ref. [2^, where a 
GaAs QD having Rq = 90 A was considered, and the 
quadratic term in Bq was neglected since only small fields 
were considered; also, only N — 2 and N = 3 occupa- 
tions were calculated, since the states were exactly built 
(not only a single Slater determinant), and the electron- 
electron interaction was included by using perturbation 
theory, justified at such radius. At zero field the energies 
for N = 2 are 16.5 meV (LS) and 16.1 meV (UHF), while 
for N = 3 they are 34.8 meV (LS) and 33.9 meV (UHF); 
so, the formalism here used indeed give smaller ground 
state energies than the LS perturbation scheme. We have 
also checked the validity of neglecting the quadratic term 
in Bq for fields smaller than 2 T. One should emphasize 
that a disadvantage of the UHF approach is that, in prin- 
ciple, it is not sure that one gets trustable information 
about the L and S expected values of QD states; on the 
other hand, the applicability of the LS scheme is highly 
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decreased as the QD occupation increases. 

We have shown how the unrestricted Pople-Nesbet ap- 
proach applied to a spherical QD under a magnetic field 
yields a reasonable description of its energetic spectrum, 
where a maximum occupation of 40 electrons has been 
considered. We have seen how both QD chemical poten- 
tial and charging energy reproduce the filling and half- 
filling structures of its energy shells at zero field. With 



the total spin expected value for each occupation in a 
given radius, we have seen that the Hund rule is satis- 
fied at zero field; however, under a finite field, we have 
shown that it is violated and, at given values of the field 
which depend on QD parameters, transitions that change 
a given ground state symmetry are observed. 
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